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ABSTRACT 

We  prove  recurrence  relations  for  a  general  class  of  multivariate  B- 
splines,  obtained  as  'projections1  of  convex  polyhedra.  Our  results  are 
simple  consequences  of  Stokes*  theorem  and  include,  as  special  cases,  the 
recurrence  relations  for  the  standard  multivariate  simplicial  B-spline. 
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SIGNIFICANCE  AND  EXPLANATION 


Because  of  their  local  support ,  finite  elements  play  an  important  role  as 
basis  functions  for  spaces  of  smooth  piecewise  polynomials.  We  have  found 
that  some  standard  finite  elements  can  be  obtained  as  •projections'  of  simple 
convex  polyhedra.  This  leads  in  a  simple  way  to  recurrence  relations  for  the 
efficient  evaluation  of  such  finite  elements. 

Even  in  the  previously  known  special  case  of  simplicial  B-splines, 
studied  in  much  detail  by  W.  Dahmen  and  C*  A.  Micchelli,  the  argument  of  the 
report  leads  to  simplifications. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


RECURRENCE  RELATIONS  FOR  MULTIVARIATE  B-SPLINES 
Carl  de  Boor1  and  Klaus  Httllig1'^ 


We  wish  to  point  out  what,  in  hindsight,  seems  obvious,  namely  that  the  recurrence 
relations  for  multivariate  B-splines  established  by  C.A.  Micchelli  [19]  and  reproved  in 
various  different  ways  by  W.  Dahmen  [6],  C.A.  Micchelli  [20],  K.  HOllig  [15]  and  H. 
Hakopian  [14]  (and  perhaps  others)  are  special  cases  of  more  general  and  very  simple 
recurrence  relations  which  are  a  simple  consequence  of  Stokes’  theorem. 

To  recall,  following  the  lead  of  I.J.  Schoenberg  [21],  the  multivariate  B-spline 

M(  • | xA,«*.,x  )  was  defined  in  [1]  by  the  rule 
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with  x  ,...,x  points  in  R11  and  conv{x„ , . . .  ,x  }  their  convex  hull,  with  vol.  (K) 
On  On  * 

the  k-dimensional  volume  of  the  set  K  ,  and 

P :  Rn — ♦  Rm:  «  — ♦  (*(1))^,  . 

Such  a  B-spline  is  a  nonnegative  piecewise  polynomial  function  of  degree  at  most  n-m  ,  its 

support  is  conv{Pxn , ♦ . . ,Px  }  ,  and  it  is  in  c11"”'1  as  long  as  the  "knots”  x 

On  u  n 

are  in  general  position. 

It  was  hoped  that  these  functions  could  be  made  to  play  the  same  basic  role  in  the 
analysis  and  use  of  smooth  multivariate  piecewise  polynomial  functions  that  their  much 
older  univariate  version  (introduced  by  Curry  and  Schoenberg  [4-5] )  had  assumed  in  the 
univariate  spline  theory.  These  hopes  have  already  borne  some  fruit?  see  Micchelli  [20] , 
Dahmen  [7-9],  Dahmen  and  Micchelli  [10-12],  Goodman  and  Lee  [13],  H&llig  [14].  The  first 
step  in  this  development  was  taken  by  C.A.  Micchelli  [19]  who  proved  the  following. 
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Theorem  1  (C.  A.  Micchelli). 

(1)  If  *  *  E  Xi**i  with  Z\  =  0  '  then 

DZM(-Ix0 . *n)  *  n  E  \M(  • 

(ii)  If  x  -  £  ^tP»i  with  £Xi  =  1  ,  then 

(n-m)  M(X|*0 . *n >  ”  "  E  . 

Here,  D^f  :■  Z  *(i)  D^f  ,  with  D^f  the  partial  derivative  of  f  with  respect 
to  its  i-th  argument.  Further,  the  equalities  asserted  in  the  theorem  must  in  general  be 
interpreted  in  the  sense  of  distributions.  In  this  connection,  Micchelli* s  starting  point 
was  the  observation  that 
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These  integrals  play  a  crucial  role  in  Kergin  interpolation  [17-19]  .  They  also  appear  in 
the  Hermite-Genocchi  formula  for  the  n-th  divided  difference. 

Consider  now,  more  generally,  a  polyhedral  convex  body  B  in  vP  ,  whose  boundary 
3B  is  the  essentially  disjoint  union  of  finitely  many  ( n-1 )-dimensional  convex  bodies  Bi 
with  corresponding  outward  normal  .  Let  M  and  denote  the  corresponding 

distributions  on  Ilf11  defined  by  the  rule 

M  $  :=  /  4>oP 

B  ,  all  test  functions  4>  • 

4>  :=  /  <J>o p 

B, 

1 

Here  ,  f  denotes  the  k-dimensional  integral  over  the  convex  set  K  in  case  K  spans  a 
k-dimensional  flat. 

Theorem  2* 

(i)  D_  M  *  -  £  <s|n.>  M.  ,  all  x  e  R° 

- -  p*  l  i  - 

(il)  (n-m)  M(Px)  =  1  <^-*111  >  M^Px)  ,  all  x  e  R°  . 
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Here,  stands  for  an  arbitrary  point  in  the  flat  spanned  by  ,  hence  the 

coefficient  <b^-x|a^>  is  simply  the  signed  distance  of  s  from  that  flat. 

The  proof  of  (i)  is  immediate: 


K,  ")<►  ”  -  /(d„,*)op  *  -  f  «>(♦•)  -  -  f  <*!«>  ♦op  • 
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As  to  (ii),  we  follow  Hakopian  [14]  who  derives  Theorem  l.(ii)  from  the  following 
B-spline  identity: 

(d-D  )m(  )  -  (n-m)M( -I*  , ) 

v  u  n  u  n 

"  . *i-r“i+r”"“n)  ' 

Here,  D  stands  for  the  differential  operator  given  by  the  rule 

k 

(Df)(x)  :»  Z  x(  j)  (D  ff  )(x) 

1-1  ’ 


for  a  function  f  of  k  variables. 
Correspondingly,  we  prove 


( iii)  DM  =■  (n-m)M  -  l  <bilili>M1 


as  follows: 
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and  this  proves  (iii)  sine*.*  <*|n^>  is  constant  on  . 
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Now,  to  prove  (ii),  conclude  from  (i)  and  (lii)  that,  for  any  z  with  P* 


x 


0  =  0  - 

®  (n-m)M(x)  -  E  <b  In^M^x)  +  Z  rzjn^M^x)  . 

Remarks,  (a)  The  convexity  assumption  is  sufficient  for  the  intended  application  but 
could,  of  course,  be  relaxed. 

(b)  Repeated  application  of  Theorem  2.(i)  shows  that  M  is  a  piecewise  polynomial  of 

degree  at  most  n-m  ,  with  possible  discontinuities  only  across  convex  sets  of  dimension 
m-1  of  the  form  P[F]  ,  with  F  a  face  of  B  .  Precisely,  M  e  with  d  the 

greatest  integer  with  the  property  that  a  d-dimensional  face  of  B  is  projected  by  P 
into  an  (m-1 ) -dimensional  set. 

(c)  This  study  was  motivated  by  the  realization  that  many  standard  finite  elements 
could  be  obtained  as  such  'projections'  of  simple  geometric  bodies  and  by  the  hope  that,  by 
using  bodies  other  than  simplices,  the  resulting  piecewise  polynomial  functions  M  might 
be  simpler  and  conform  more  easily  to  standard  meshes.  First  results  along  these  lines  are 
contained  in  121  and  13]. 
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